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Abstract. 

We show that if the primordial classical perturbations were generated by the gravitational 
particle creation during inflation, and followed by an evolution of quantum-to-classical tran- 
sition, the time dependent correlation of these perturbations is long-tailed with a correlation 
time larger than the Hubble-time. Consequently, the inflationary perturbations are locally 
scale-scale correlated. Hence, the interaction of the fields during inflation can be explored via 
the detection of the local scale-scale correlation of the CMB fluctuations. 



It is fundamentally important to understand the nature of the primordial cosmological 
perturbations which are responsible for the formation of structures that we see today. The 
inflation paradigm for the early universe assumes that the cosmic perturbations are initiated 
by the process of particle creation from vacuum in the background gravitational field of the 
expanding universe |l]] . The subsequent decoherence of the quantum fluctuations leads to a 
quantum-to-classical transition and gives rise to the initial perturbations in the radia- 
tion dominated era on scales beyond the Hubble radius. Therefore, the earliest evolution of 
the inflationary cosmological fluctuations is governed by the gravitational particle creation 
accompanying the quantum-to-classical transition. Recently, the mechanism of the gravita- 
tional particle creation is re-emphasized with the development of the quintessential inflation 
model However, it still seems to lack some testable predictions of these mechanisms. 

In this Letter, we will show that the time-dependent correlation of the primordial perturba- 
tions is one of such testable predictions. In the "standard" inflation, i.e. the slow-roll inflation 
caused by a single scalar inflaton <j>, the correlation time of the primordial perturbations is 
longer than the Hubble time H . This correlation will entail a non-trivial observable feature 
- the local scale-scale correlation of the initial perturbations. 

The time-dependence of cosmological perturbations is generally not observable in non- 
inflationary models, because the observed mass field does not contain two or multi-time in- 
formation of the perturbations. In contrast, the inflationary scenario provides an excellent 
prospect to observing the time-dependence of the initial perturbations in a given region on 
the scale of horizon. Owing to the so-called "first out - last in" feature, a fluctuation crossing 
over the Hubble radius at the moment t x will yield a perturbation at the end of the 

inflation, t / , with the physical wavenumber given by Q 

kpoce 11 ^-^. (1) 
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Consequently, the longer the physical length of the perturbations, the earlier the time of be- 
coming the super-horizon scaled, and vice versa. Therefore, if perturbations with the same 
physical scale at two different instants t* and t% are correlated during the super-horizon evo- 
lution, the perturbations of the corresponding physical scales fcp(i* ) and kp^) at the same 
time tf will be correlated. Equation (1) prescribes exactly the mapping of the correlation be- 
tween two perturbations on the same physical scale with different horizon-crossing moments 
to that of two perturbations at the same time with different physical scales. The latter is 
actually observable. If the length scales of perturbations are larger than the Hubble radius 
at the decoupling era, the perturbations were less contaminated by the post inflationary 
reheating and other small scale processes. Thus, one can expect that the two-scale correla- 
tion of the cosmic background temperature fluctuations retain the information of the time 
dependent behavior of the initial perturbations. 

The time-dependent correlation is not a new issue. It is generally believed that the per- 
turbations originated from the quantum fluctuations of vacua during inflation have large cor- 
relation time since they vary slowly in the super-horizon regime. However, the rough idea of 
the lengthy correlation time is inadequate to establish the detectability of the time-dependent 
correlation , and one must study the temporal correlation of the perturbations in both quan- 
tum and classical regime to reach the equation suitable of calculating the time correlation in 
practical. 

Let us calculate the time-dependent behavior for the "standard" slow-roll inflation gov- 
erned by a real massless scalar field <f> described by the action 



where V(<f>) is a self-interaction potential. Under the slow-roll condition, the time-scale of the 
self-interaction is much longer than 1/H. That is, the self-interaction is negligible if we study 
only the evolution on time scales smaller than the duration of inflation. Without interaction, 
the evolution of perturbations of the free cj) field mode will be coherent. Thus, the <f> field 
self-interaction is ineffective of violating the temporal coherence during inflation. 

For the simple model Eq. (2), the dominant interaction during inflation is the coupling 
between the <fi field and the gravitation of the expanding universe, governed by the gravita- 
tional scalar density \J—g in the integral. Therefore, in the "standard" inflation model, the 
origin and evolution of the temporal coherence of the primordial perturbations are actually 
determined by the dynamics of the gravitationally driven particle creation. 

In a de Sitter space, the scalar field cf>(r, t) can be described as a superposition of free modes 
with coordinate and conjugate momentum variables given by 0(k) = (2tt)~ 3 / 2 J 0(r)e _lk r c? 3 r 
and 7r(k) = a 2 <j>'(k.), respectively, where k is the comoving wavevector and a — e Ht represents 
the cosmic scale factor || . The notation / denotes the derivative with respect to the conformal 
time r = J dt/a(t) = —(l/H)e~ Ht . Therefore, the time-dependent behavior of the <fi field can 
be studied via the modes k. 

The quantum nature of the 4> field is specified by the equal-time commutation relation 
given by [0(k, r), 7T (k, t)] — i<5^ 3 ^ (k — k'), which yields 




(2) 



0(k,r) 
?r(k, r) 



CT k (r)a k + Ck(r)aL k , 
a 2 cr k (r)ak + a2 °k'( T ) a t k - 



(3) 



The a k and are respectively, the annihilation and creation operators, satisfying commuta- 
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tion relations [ok, a^,] = 5^ 3 ^(k — k'). The time dependence of the mode k is given by 

Ht-i^- e' lkT , (4) 



o"k(r) 



'2k 

If the system is in the vacuum state at a given time To, we have 

a k |0,r )=0 (5) 

for all k. In the Schrodinger representation, the evolution of the system leads the vacuum 
state at tq developing into |0,t) = S\0,To) at a later time r, where S is the S'-matrix. The 
state |0, r) is no longer a vacuum state. The average number of the particles gravitationally 
created in mode k is N^t) = (0, To|S , ~ 1 aj c a_k>S'|0, r ). 

One cannot directly calculate the time-dependent correlation of classical density pertur- 
bations by the state |0, r), because at a different r the states are quantum coherent. We 
should first find the condition of quantum decoherence for the system at different times. In 
the coordinate representation 0(k), the Schrodinger wave function of the state |0,t) is given 
by the product of the wave functions for all k modes : 

*[0(k), r] = 1 exp {- JML[1 - l 2F(k, r)]\ , (6) 
vV|CT k (r)| L 2|cr k (r)| 2 J 

where 

Equation (6) is actually the ground state wave function of a harmonic oscillator with coordi- 
nate |0(k)| which possesses a time-dependent variance in itself. 

Thus, the quantum coherence between the coordinates at different times can be calculated 

as 

(*(r 2 )|^(k)^(k)|*(r 1 )) = \ (8) 



; [l-i2F(k,n)] + - ) .,„ [l + i2F(k,75)] 



2| ( r k (r 1 )| 2 L v ' 1/J ' 2|a k (r 2 )| 2 



-3/2 



When i ^> i? -1 or \Ht\ = e~ Ht 1, the term Ht in Eq. (4) can be ignored, and |tTk(T)| 
becomes r-independcnt. Accordingly, |CTk(7"i)| — Ick^)!) an d the ratio of the two-time matrix 
element to the diagonal element reduces to 



\ e H tl _ e m 2 



(9) 



Evidently, Eq. (9) is quickly approaching zero which implies that the decoherence of the 
system for all \ti — 1 2 | > l/H, once ke~ Htl and ke~ Ht2 are less than H . Thus, the relationship 
between the k-mode fluctuations at times t\ and t 2 can be treated as classical perturbations 
provided that both physical wavelengths, e Htl jk and e Ht2 /k, are super- Hubble sized with the 
time separation \ti — i 2 | being greater than one Hubble time. 

The quantum decoherence does not exclude the classical time correlation of perturbations 
M . The classical correlation can be derived by means of the Wigner distribution functions of 
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the system. For the k mode, it reads 

W(<f>(k), 7r(k), r) = / d (£-) e"^* (0( k ) - |, r) * (^(k) 



= iV k ^exp 



K(r)| 



exp 



-kk(r)| 2 



TT(k)- 



f(k,r) 
kk(r)| 2 



«£(k) 



(.10) 



Since |crk(7-) | is r-independent when t ^> 1/H, the first exponent on the right hand side of 
Eq.(10) represents the Gaussian probability distribution of the fluctuation of the scalar field 

The time-dependent behavior of the classical perturbations can be extracted through the 
second Gaussian probability distribution in Eq. (10). In the classical limit Ti — > 0, it yields 



W(0(k),7T(k),' 



iVk exp 



\m\ 2 

kk(r)P 



^ (k) -RW 0(k) 



This shows that the classical trajectory of the k mode in phase space is given by 

F(k,r) 



TT(k) 



kk(r) 



(k) = 0. 



(11) 



(12) 



The degree, or the effectiveness, of classical correlation can be measured by the relative sharp- 
ness of the classical trajectory in phase space, which is defined as the ratio of the dispersion 
in momentum to the magnitude of the average of the momentum 1 7r(k) | ||. By virtue of Eq. 
(10), this ratio is 



l/kk(r)| 



1 



KF(k,r)/|<r k (r)| 2 )^(k)| 



FCk,r) 



«1, if \kr\ = (k/H)exp(-Ht) -C 1. (13) 



Therefore, the classical perturbations of super-Hubble size, i.e. keT m <C H, are perfectly 
coherent. Namely, after the quantum-to-classical transition, the evolution of the initial per- 
turbations can be calculated by the classical trajectory Eq. (12). 

The scale t c of the time-dependent correlation between classical perturbations of k mode 
can be characterized by 



1 



d<f>(k)/dt dr 7r(k) 



t c 0(k) 9ta 2 0(k) 

Hence, for super-Hubble sized perturbations, we have 

t c H > 1. 



k 2 e -2Ht 

H ' 



(14) 



(15) 



That is, the correlation time of the primeval perturbations is much longer than the Hubble 
time H- 1 . 

This result is anticipated because the only interaction of the inflation in the "standard" 
model during inflation is the gravitational particle creation. After the quantum-to-classical 
transition, the perturbed scalar field consists of classical free motion. Since there is no cou- 
pling between either comoving modes k or physical modes k p , the evolution of the classical 
perturbations of the <fi field are coherent. Therefore, one may conclude that the large cor- 
relation time of the initial perturbations probably is a generic feature for models with free 
inflaton, i.e., besides the gravitational coupling of the inflaton, there are no interactions and 
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self-interaction during the slow-roll phase. This point can be seen more clearly if we establish 
an equation for the classical perturbations of the <fi field. 

Actually, one can regard Eq.(14) as an evolution equation of the ensemble averaged </>(k), 
i.e. d(j}(k)/dt ~ — (fc 2 e~ 2H */iJ)</)(k). The uncertainty of this equation can be estimated by 
the time scale given by 



1 



Ad(f>(k)/dt 



|A7r(k)| ^ k 3 e~ 3Ht 



We have t' c 3> t c and therefore, the classical evolution equation for the k mode perturbations 
in phase space is reasonable in average. 

The two equations (14) and (16) of the ensemble averaged mode <fi(k) can be combined 
into the following Langevin-like equation for stochastic variable <j>(k) 

where the small Gaussian noise r/^ has zero mean with variance ~ k 3 / 2 e~ 3Ht /H. Therefore, 
the ensemble average of (17) yields Eq.(14), and the variance gives (16). Since Eq.(17) is linear 
to (f>(k, t) , it is straightforward to get the equation of <p(r, t) in the comoving r- representation. 
After employing the physical coordinates in which dxi — e Ht dri (i = 1, 2, 3), Eq.(17) becomes 

*^~iv^(x,t) +,,(«, t), (18) 

where V 2 stands for the Laplacian of the physical coordinates x. Because the variance of 
the stochastic noise is scaled as (|r?(x, t)\ 2 ) oc e~ 3Ht , the noise term is strongly suppressed as 
Ht > 1. 

Equation (18) shows that the correlation time for a mode with the physical wavevector k p , 
i.e. </)(x, t) ~ e _lkp X , is proportional to kz 2 . This is typical to the diffusion-like soft modes in 
the theory for non-equilibrium systems || . If a perturbed field consists of the superposition of 
such soft long wavelength perturbations described by Eq.(18), the time dependent correlation 
function of the perturbed <f> field possesses a long-tail [[[o) . This behavior is more transparent 
when considering modes characterized by e - lwt -* k P x . For these modes Eq.(18) gives rise to 
the dispersion relation as ui = —ik 2 /H. Apparently, the relaxation time — —> oo when 

k p — > 0, and these modes lead to long-range correlations in general. 

Since we are only interested in the perturbations consisting of modes with physical wave- 
length larger than if -1 , the field ^(x, t) crossing over the horizon can be prescribed by 

<P(x,t) = f d 3 k p $ kp e- 1 ^-^, (19) 

Jk p <H 

where <f>\^ p is the amplitude of the perturbation for modes k p . Thus, the time-dependent 
correlation function is determined by 

(0(x,ii)0*(x,t2))= jf ^d 3 k p (0 kp ^)exp (^-h\ tl -t 2 \\ . (20) 

At a given instant, the distribution of (4>4>*) with respect to the k p modes should be the same 
as that along the k modes, since the relation between k p and k is specified by a constant 
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multiplier. On the other hand, Eqs.(4) and (11) imply that, (|0(k)| 2 ) ~ |ok(T)| 2 cx k 3 in the 
super-horizon region. Therefore, (</>k p 0k ) oc fc~ 3 , and the integral (20) is infrared divergent. 
However, the duration of the inflation spans only a finite period of time, there must be an 
infrared cutoff at k ~ H/R with R^$> 1. Thus, (<pk p 4>^ ) can be approximated as a quantity 
independent of k p , and Eq.(20) yields 



(^(x,fi)^(x,ta))« I | ti _ ta | j > if |*i >*o, (21) 

where the correlation time t c ~ C/H, and the constant C > 1, i.e. the correlation time 
is longer than the Hubble time. Hence, during the super-Hubble evolution, the field <j> is 
coherent regardless of the horizon-crossing moments of the perturbations provided that time 
difference is less than C/H. This result illuminates that the gravitational particle creation can 
be considered as a dissipation source for the non-equilibrium inflaton <p during the slow-roll. 
The relaxation of that dissipative process is dominated by soft modes and subsequently gives 
rise to the classical perturbations with a long-tailed time correlation. On the other hand, the 
soft modes do not induce any long range spatial correlation and therefore, the perturbations in 
different spaces x retain themselves. That is, the random field of the inflationary perturbations 
is "double-faced" : its spatial distribution is Gaussian, while the temporal or scaled distribution 
is coherent. This property leads to a measurable effect if we use field variables based on a 
space-scale decomposition. 

Let us consider such space-scale decomposition given by the bases x (x') which are 
localized, orthogonal and complete. The indices k p and x denote a cell in the phase (x, k p ) 
space, i.e. k p — > k p + Ak p , x — > x + Ax, and the differential volume element AxAk p ~ 27r. 
Wavelet analysis provides various bases for this sort of space-scale decomposition jl]J . With 
a suitable wavelet decomposition, the (f> field can be described by the variables 6<pk p ,x defined 
as 

(x')* V -( x ')dx'. (22) 

Obviously, S4>k p ,x represents the amplitude of the perturbation of mode (k p ,x), i.e. the fluc- 
tuations of the field at position around x and on scale around k p . 

According to Eq.(l), the amplitude 54>k p ,x is mainly attributed to 0(x) with the horizon- 
crossing time corresponding to k p . Consequently, the time-time correlation of Eq.(21) implies 

{SKuMk**) ^ 0, (23) 

if the difference between the horizon-crossing times corresponding to k p \ and k p i is less than 
C/H. Equation (23) is a local (at spatial area around x) scale-scale {k p \ and k P 2) correlation. 

By means of the variables &4>k p ,x one can easily perceive this "double- faces" feature of the 
inflationary perturbations. The variables S(j)k p .x which are Gaussian with respect to x entail 

{5<pk p ,x8<j)k p ,x>) = P(kp)5 XjXf , (24) 

and all higher order cumulants of §4>k p ,x are zero, in which P{k p ) is the power spectrum of the 
Gaussian fluctuations 5tfik p ,x jjjj- Meanwhile, the variables at different scales k p may have 
carried some correlation. As a simple example, let us consider 

8<t>k' p ,x = ot5cj) kpiX , (25) 

where a is a constant. In this case, the perturbation 5<f>k' ,x with respect to x is also Gaussian, 
i.e. its variance is P{k' ) — a 2 P(k p ) and all higher order cumulants of 5(f>k< , x are zero. However, 
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the scale-scale (k'-kp) correlation at the place is significant, i.e. (6(j>k> ,x^k p ,x) ^ ol. Therefore, 
a Gaussian power spectrum P(k p ) can coexist with a local scale-scale correlation [Eq.(23)]. 
It should be emphasized that as a statistical measure, the local scale-scale correlation is 
independent of the Gaussian power spectrum. One cannot determine whether there is and/or 
how strong it is the local scale-scale correlation directly by their power spectrum. 

Recently, the search for the local scale-scale correlation of large scale structure samples has 
attracted attention. Using the wavelet technique, the local scale-scale correlations of QSO Lya 



forests 1 13 and COBE data jlj] have been studied. On the other hand, the time-dependent 
correlation of the primordial perturbations is sensitive to the interaction of the inflationary 
field. It is also sensitive to dissipation during the inflation. For instance, if the cosmic inflation 
is thermally dissipated |l5| ], the correlation time of the initial cosmic perturbations will be 
significantly changed by a thermal-dissipative term in Eq.(17). Hence, the local scale-scale 
correlation behavior of the CMB fluctuations and other relevant samples would be useful to 
discriminate among models of inflation. 

* * * 
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